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Abstract 



Numerical simulations of conduction through a disordered microbridge be- 
tween a normal metal and a superconductor have revealed an anomalous in- 
sensitivity of the conductance fluctuations to a magnetic field. A theory for 
the anomaly is presented: Both an exact analytical calculation (using random- 
matrix theory) and a qualitative symmetry argument (involving the exchange 
of time-reversal for reflection symmetry). 
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Universal conductance fluctuations (UCF) are a fundamental manifestation of phase- 
coherent transport in disordered metals |l]J^. The adjective "universal" describes two aspects 
of the sample-to-sample fluctuations of the conductance: (1) The variance varG is of order 
(e^//i)^, independent of sample size or disorder strength; (2) varG decreases precisely by a 
factor of two if time-reversal sjTumetry (T) is broken by a magnetic fleld. The universality of 
this factor of two has been established both by diagrammatic perturbation theory ||l|,0] and 
by random-matrix theory In the former approach, one has two classes of diagrams, 

cooperons and diffusons, which contribute equally to var G in the presence of T. A magnetic 
fleld suppresses the cooperons but leaves the diffusons unaffected, hence varG is reduced 
by |. In the latter approach, the universality of the factor-of-two reduction follows from 
the Dyson- Mehta theorem 0, which applies to the variance var A of any observable A = 

a{Tn) which is a linear statistic on the transmission eigenvalues Tn [||. The crossover from 
a linear to a quadratic eigenvalue repulsion upon breaking T directly leads to a reduction 
by I of var A p. 

The situation is qualitatively different if the normal-metal conductor (N) is attached 
at one end to a superconductor (S). At the NS interface the dissipative normal current is 
converted into a dissipationless supercurrent, via the scattering process of Andreev reflection 
T0| : An electron incident from N is reflected as a hole, with the addition of a Cooper 



pair to the superconducting condensate. The conversion from normal to supercurrent has 
essentially no effect on the average conductance, provided that the interface resistance is 
negligibly small ||rT|. However, the effect on the conductance fluctuations is striking: The 
variance is still universally of order (e^//i)^, but it has become insensitive to the breaking 



of T. Numerical simulations by Marmorkos, Jalabert, and one of the authors |T2[ of a 
disordered wire attached to a superconductor have shown that the variance is unaffected 
by a T-breaking magnetic fleld, within the 10% statistical uncertainty of the simulations. 
This does not contradict the Dyson-Mehta theorem, because the conductance Gns of the 
NS junction is a linear statistic in the presence — but not in the absence of T |]T3[. One 



wonders whether there is some hidden symmetry principle which would constrain varCws 
to be the same, regardless of whether T is broken or not. No such symmetry principle has 
been found, and in fact we do not know of any successful generalization so far of the theory 



of UCF to quantities which are not linear statistics [|14 



Here we wish to announce that we have succeeded in the analytical calculation of var G^s 
in the absence of T, using techniques from random-matrix theory. We flnd that varGNs 
for a disordered wire attached to a superconductor is reduced by (2 — OO/vr^)"^ ^ 0.929 
upon breaking T. This number is sufficiently close to 1 to be consistent with the numerical 
simulations and sufficiently different from 1 to explain why attempts to flnd a rigorous 
symmetry principle had failed. Still, we have been able to flnd an approximate symmetry 
argument, which explains in an intuitively appealing way why the number (2 — 90/7r^)~^ is 
close to 1. Our theory is more generally applicable than to a disordered wire: It applies to 
any NS junction for which the probability distribution P{S) of the scattering matrix S of 
the normal region depends only on the transmission eigenvalues T„. (Such a distribution is 
called "isotropic" 0.) As two examples we consider a disordered metal grain and a ballistic 
constriction in a disordered wire. 

Starting point of our calculation is the general relation between the conductance of the 
NS junction and the scattering matrix S of the normal region ||T3], 
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GNs = 2Gotrmm^ Go = 2e^/h, (la) 
m = Vf{l + uVRu*^)-^uVf, u = w^wl (lb) 

We used the polar decomposition 

^_(vi \ (iVR VT \(v2 \ 

where vi, V2, wi, and W2 are N x N unitary matrices (A^ being the number of propagating 
modes at the Fermi level in each of the two leads attached to the sample). The matrix T 
is a diagonal matrix with the N transmission eigenvalues Tj G [0, 1] on the diagonal, and 
R = 1 — T. In the presence of T, one has 5* = S^, hence W2 = wj, hence u = 1. (The 



superscript T denotes the transpose of the matrix.) Eq. (1) then simplifies to [T^ 

GNs(r)=2GoEnT„'(2-T„)-2, (3) 

and var^Ns follows directly from general formulas for the variance of a linear statistic on 
the transmission eigenvalues P,p^. In the absence of T no such simplification occurs. 

To compute var^Ns = (G^ns) ~ (G^ns)^ in the absence of T, we assume an isotropic dis- 
tribution of S, which implies that the average (• • •) over the ensemble of scattering matrices 
can be performed in two steps: (■■■) = ((• • ■)u)t, where (• ■ ■)„ and (■ ■ ■)t are, respectively, 
the average over the unitary matrices u and over the transmission eigenvalues Tj. It is con- 
venient to add and subtract ((GNs)n)r; ^o that the variance of the conductance splits up 
into two parts, 

varGNs = {{G^s)1)t " {{G^s)u)l + 

{{Gls)u-{G^s)l)^, (4) 

which we evaluate separately. 

The first part is the variance of (Gns)u over the distribution of transmission eigenvalues. 
As a consequence of the isotropy assumption, the matrix u is uniformly distributed in the 
group U{N) of N X N unitary matrices [Q. To evaluate (G'ns)u we need to perform an 
integral over U{N) of a rational function of u, according to Eq. (1). Such matrix integrals 



are notoriously difficult to evaluate in closed form []TB[, but fortunately we only need the 



large- limit. Creutz and Mello [|T8| have given general rules for the integral over IA{N) 



of polynomial functions of u. By applying these rules we find that 

{ti T{uVRu*^YuTu\^u^VRu^f)u = 
6„NT',{l-nfP + 0{l), (5) 

where we have defined the trace = A^~^ UtT^ . It follows that, up to corrections of order 
unity, 

(Gns)« = 2G,N±tI{1 - nf^ = (6) 
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Since Tk is a linear statistic, we know that its fluctuations are an order smaller than 
the average [^. This implies that, to leading order in var/(rfc) = [/'(r^)]^ var r^.. The 
variance of Eq. (^) is therefore 

Note that the leading term in Eq. (0) is 0{1). 

We now turn to the second part of Eq. (^), which involves the variance (G^g)^ — (^ns)u 
of Gns over U{N) at flxed transmission eigenvalues and subsequently an average over the 
Tj's. The calculation is similar in principle to that described in the preceding paragraph, 
but many more terms contribute to leading order in Here we only give the result, 

{{Gls)u - (Gns)')^ = lec^ (2 - in))-' {A{nf - 8{nf + 9{n)' - 4{n)' + 2{n)'- 

^r,)\r2f - 4(n)(r3) + ^{r,)\r-^ - 2{r,f{r,)] + 0(l/iV). (8) 

The sum of Eqs. and (H) equals varGNS; according to Eq. (p. The resulting expression 
contains only moments of the transmission eigenvalues. This solves the problem of the 
computation of var Gns in the absence of T, since these moments are known. 

For the application to a disordered wire (length L, mean free path one has the variance 
@J|] A^Varri = ^, and averages [|l9l (r^) = \{i j L)Y {k) jV {k ^ |). Substitution into 



Eqs. (0) and (|) yields (in the diffusive limit IjL 0) 

var Gns (no T) = ^G^ ^ 0.533 G^. (9) 

This is to be compared with the known result in the presence of T |T^ 

var GNs(r) = (if - 487r-^)G2 ^ 0.574 G^. (10) 

Breaking T reduces the variance by less than 10%, as advertised. 

We would like to obtain a more direct understanding of why the two numbers in Eqs. (^ 
and ([T0|) are so close. To that end we return to the general expression (1) for the conductance 
Gns of an NS junction, in terms of the scattering matrix S of the normal region. We compare 
Gns with the conductance Gnn of an entirely normal metal consisting of two segments in 
series (see Fig. |^). The flrst segment has scattering matrix S", the second segment is the 
mirror image of the flrst. That is to say, the disorder potential is specularly reflected and 
the sign of the magnetic fleld is reversed. The system NN thus has a reflection symmetry 
(iS), both in the presence and absence of T. The scattering matrix of the second segment 
is ES'S, where S is a 2N x 2A^ matrix with zero elements, except for Sj^7v+i = "^nma — 1 
(z = 1, 2, . . . , A^). (The matrix S interchanges scattering states incident from left and right.) 
The conductance Gnn follows from the transmission matrix through the two segments in 
series by means of the Landauer formula, 

Gnn(5) = GotrmW^ (11a) 
m! = Vt{1 + u'\/Ru'Vr)-^u'Vt, u' = W2W1. (lib) 



4 



The difference between Eqs. (1) and (|T^) is crucial in the presence of T, when W2 = wf, so 
that u = 1 while u' is some random (symmetric) unitary matrix. However, in the absence 
of T, wi and W2 are independent, so that both u and u' are randomly distributed unitary 
matrices. We have repeated the calculation of the variance starting from Eq. (0), and found 
that vartr mm^ = vartr m'm'^, hence 

varGNslno T) = 4varG'NN('5, no T). (12) 

The system NN is special because it possesses a reflection symmetry. Breaking S amounts 
to the replacement of the mirror-imaged segment by a different segment, with scattering 
matrix S' which is independent of S but drawn from the same ensemble. Breaking S reduces 
the variance of the conductance fluctuations by a factor of two, regardless of whether T is 
present or not, 

varGNN('5) = 2varG'NN(no S). (13) 

We have checked this relation by an explicit calculation, but it seems intuitively obvious 
if one considers that the eigenstates separate into even and odd states which fluctuate 
independently. Since breaking T by itself reduces the variance of Gnn by a factor of two, 
we may write 

varG'NN('5, no T) = varG'NN(^, no S). (14) 

Eqs. (p!^)-([T^) are exact, and hold for any isotropic distribution of the scattering matrix. We 
need one more relationship, which is approximate and holds only for the case of a disordered 



varG'Ns(^) ~ 4varGNN(T, no S). (15) 

Taken together, Eqs. (p^-(p^ imply the approximate relationship varGNs(^) ~ 
varGNs(iio T). The exact calculation shows that the approximation is accurate within 10%. 
We now understand the insensitivity of the conductance fluctuations to a magnetic field as 
an exchange of symmetries: Breaking T induces S, thereby compensating the reduction of 
var Gns- 

We have emphasized the general applicability of Eqs. (0), (||) and (p!2D-(p!^), which hold 
not just for a disordered wire, but for any isotropic distribution of the scattering matrix. 
We illustrate this by two examples. The first is an NS junction consisting of a disordered 
metal grain (see Fig. ^ja). The coupling of N and S to the grain occurs via ballistic point 

we 



contacts (width much smaller than the mean free path in the grain). Following Ref. pl| , 
may assume that the scattering matrix of the grain is distributed according to the circular 
ensemble of random-matrix theory. This is an isotropic distribution. The relevant moments 
of the transmission eigenvalues in the absence of T are [^] (r^) = 7r~^/^r(fc + ^)/T{k + 1), 
A^^varri = ^. Substitution into the general formulas (0) and (^) yields 

varGNs(no T) = iffG^ ^ 0.5276"^, (16) 



which is again close to the known result in the presence of T E^, 
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var Gm{T) = ^Gl ^ 0.563 Gl (17) 

The second example is a ballistic constriction (point contact) in a wire which is connected 
to a superconductor (see Fig. ^). The point contact has conductance NqGq, which we 
assume to be much smaller than the conductance Ni/L of the disordered wire by itself. As 
discussed in Ref. we may assume an isotropic distribution of the scattering matrix of the 
combined system (point contact plus disordered wire). The moments of the transmission 
eigenvalues are ^ (rfc) = Nq/N, N'^vht n = 0{NoL/Ni)'^. Substitution into Eqs. (0) and 
(I) yields, in the hmit NqL/NI 0, 

varG'Ns(no r) = iG^. (18) 
In contrast, if T is not broken, the conductance fluctuations are suppressed in this limit 



22,23 : 



varGNs(T) = OiNoL/Nif < Gl (19) 

In this geometry a magnetic field greatly enhances the conductance fluctuations. The reason 
that a disordered wire with a constriction behaves so differently from an unconstricted wire, 
is that the relation (|1^) does not hold in the presence of a constriction. However, the 
general relationship (|I2|) does hold, and indeed the result (0) is four times the variance of 
a structure consisting of two point contacts in series with a reflection symmetry. 

In summary, we have solved the problem of universal conductance fluctuations in normal- 
metal-superconductor junctions in a magnetic fleld, under the assumption of an isotropic 
distribution of the scattering matrix of the normal region. We flnd that breaking time- 
reversal symmetry introduces an approximate reflection symmetry into the structure of the 
scattering matrix. This reflection symmetry compensates the reduction of the conductance 
fluctuations due to breaking of time-reversal symmetry, and explains the anomalous insen- 



sitivity of the fluctuations in a magnetic fleld discovered in computer simulations [|T2|. We 
know of no other situation in physics where application of a magnetic fleld creates rather 
than destroys a spatial symmetry. 

Discussions on this problem with A. Altland are gratefully acknowledged. This work was 
supported by the Dutch Science Foundation NWO/FOM. 
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FIGURES 



FIG. 1. (a): Schematic drawing of a disordered normal metal (N) connected to a supercon- 
ductor (S), in a time-reversal symmetry (T) breaking magnetic field B. In (b) the normal region 
is connected in series with its mirror image. As indicated, the magnetic field B changes sign upon 
refiection. The variance of the conductance fluctuations in (a) is exactly four times the variance 
in (b). The variance in (b) is exactly two times the variance in the absence of the reflection sym- 
metry {S). The exchange of T for S explains the insensitivity of the conductance fluctuations to 
a magnetic field, as discussed in the text. 

FIG. 2. (a): Schematic drawing of an NS junction consisting of a disordered metal grain 
(shaded), (b): A disordered normal-metal wire (shaded) containing a point contact. 
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Figure 1 . 




